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ABSTRACT. A paradigm model is suggested for describing the 
diffusive limit of trajectories of two Lorentz disks moving in a 
finite horizon periodic configuration of smooth, strictly convex 
scatterers and interacting with each other via elastic collisions. 
For this model the diffusive limit of the two trajectories is a mix- 
ture of joint Gaussian laws (analogous behavior is expected for 
the mechanical model of two Lorentz disks). 
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1. Introduction 

Beside the dynamics itself, the joint motion of two particles in- 
teracting with each other and with a dynamical environment also 
depends on the spatial dimension. The first model where this ques- 
tion was addressed (cf. [Sz 80]) was a one-dimensional mechanical 
one. There, the asymptotically diffusive motions of the two parti- 
cles were either glued together or were independent depending on 
the initial distance of the particles. The model was actually that of 
Harris and Spitzer, (see [S 69J) (equilibrium dynamics of elastically 
colliding point particles) generalized by Major and Szasz, [M Sz 801 
(non-equilibrium dynamics). In a related one-dimensional - random 
- collision system (cf. [KL PS 8311 ) the joint motions were dependent. 
On the other hand, Kipnis and Varadhan ([KV 86]) have shown that 
the diffusive limits of two particles in a symmetric exclusion pro- 
cess on Z d with d > 1 are independent Brownian motions except the 
one-dimensional nearest neighbor case when the motions are subd- 
iffusive. 

It is worth noting that the joint motion of two particles got also 
studied in the physics literature, e. g. the mutual dynamics of pairs 
of atoms in a dense Lennard-Jones liquid in IP BV 841 . 

Returning from stochastic dynamics to a deterministic one, let us 
consider the planar, finite-horizon Lorentz process with a periodic 
configuration of scatterers. It is known that its limit in the diffusive 
scaling is a Brownian motion (cf. |BS 81 j and IIBCS 911). Of course, 
two point like Lorentz particles do not interact, but if we take two 
small disks then the case is quite different. 

A simple fact: The motion of one small disk is itself isomorphic 
to a Lorentz process, so its diffusive limit is again the Wiener pro- 
cess. However, if one considers two small Lorentz disks, then the 
naive heuristics would suggest that, since the two particles collide 
very rarely (i. e. O(logn) times during the first n collisions), the 
situation is similar to the locally perturbed Lorentz process where 
the diffusive limit is the same Brownian motion as it was for the 
unperturbed Lorentz process (cf. [D SzV 09l ). This analogy is, how- 
ever, misleading and the aim of the present work is exactly to clarify 
the situation. The difference with the preceding models is the inter- 
action: elastic collision of the disks also changes the energies of the two 
particles. Moreover, in dimension two, by borrowing heuristics from 
random walk theory (cf. US 761 ) and estimates from [DSzV 08], one 



can convince himself /herself that the time intervals between consec- 
utive collisions have a slowly varying tail. Consequently, for large t, 
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the last collision of the disks preceding t befell at time o(t) with a 
probability close to one. Thus the energies of the disks at time t de- 
termining the limiting variance are the random energies obtained at 
the aforementioned last collision before t, ergo the diffusive limit of 
each disk is a Brownian motion with a random covariance (and their 
joint limit can already be calculated based upon the previous line of 
ideas). 

The goal of the present work is to make the above heuristic argu- 
ment precise on the level of a stochastic model mimicking the deter- 
ministic model of two Lorentz disks. 

Our model is, roughly speaking, a colliding system of two random 
walks with internal states where the speeds of the particles are rep- 
resented by exponential clocks and are included in the set of inter- 
nal states. The model, our main conditions and our main result are 
described in subsections 2.1-2.2. Subsection 2.3 contains the main, 
often new, probablilistic concepts and results and a sketch of our 
proof. Section 3 is devoted to the verification of our local limit the- 
orem for general random walks with internal states and some corol- 
laries, whereas section 4 the the proof of our main Theorem 1. Sec- 
tion 5 contains some remarks. Finally, the proofs for our results for 
Markovian renewal processes and scaled type Markovian renewal 
processes is provided in the appendices. 

2. The Model. Main result. Methods 

2.1. The Model. The dynamics of two Lorentz disks will be mod- 
eled by two continuous random walks with continuous internal states 
whose steps are independent whenever the walkers are at different 
lattice sites. If they are on the same site, then their interaction is given 
by a collision operator (see below). 

2.1.1. Continuous Time Random Walks with Continuous Internal States. 
Discrete time random walks with a finite number of internal states 



were introduced by Sinai, MS 8111 where the internal states were meant 
to represent elements of a Markov partition. The theory was elab- 
orated in a series of works HKSz 8 3 , KSz 84, IKSSz 86 ll. In our case 
the internal states will also represent particle velocity therefore we 
have to consider random walks with internal states where the inter- 
nal states belong to a more general state space. Moreover, for being 
able to include speed we take continuous time. In MKSz 831 , a lo- 
cal limit theorem was established for random walks on Z d with a 
finite number of internal states and we will also use much of the 
techniques presented there. 
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Definition 1. (Sinai, 1980) Let H, \H\ < oo be the set of states. On the set 
Z d x H the Markov chain £ M = (rj n , e n ) is a Random Walk with Internal 
States (RWwIS) if for V x n , x n+1 eZ d , u n € H,Ac H 

P(£n+1 = (x n+1 ,U n+1 ),U n+1 G A\l, n = (x n ,U n )) = p Xn+1 -x n (u n , A) 

Of course, {z n ; n > 0} n is also a Markov chain due to the spatial transla- 
tion invariance. 

Our paradigm for the mechanical model will be introduced in two 
steps. First, the individual motion of each of the two particles will be 
a continuous time Markovian random walk with internal states (we 
will abbreviate it by RWwIS again; it will always be obvious whether 
we are talking about the discrete or the continuous time case) with 
some general state space H and a constant exponential-jump rate 
A > 0. (So far we do not specify H). We just note, however, that later 
A will be included among the internal states of the full two particle 
system to permit its change at collisions of two particles. 

Definition 2. Assume we are given a rate A > and a family 

(l) {P x (v,.)\xez d \{0}} 

of substochastic kernels over H such that Q = E x ez d \{o} ^ x zs a stochastic 
kernel over H. A continuous time pure jump Markov process £f = (rjt, £f) 
with state space Z d x H - is called a (generalized) Random Walk with In- 
ternal States (RWwIS) if 

F(.£t+dt = £t) = l-Mt + o(dt) 

and for every (x t , u) G Z d x H, VA C H and x t+ at — Xt ^ 

*(£t+dt = (x t+dt ,v'),v' <e A\Ct = (x t ,v)) = \P Xt+dt - Xt (v,A)dt + o(dt) 

In other words, the kernel for a jump to x € Z d \ {0} is described 
by 

Pxf{v) = l f(v')P x (v,dv'), f e Loo(H) 

J H 

and the transition operator for the discrete time Markov chain {e n } n >o 
of subsequent internal states is Q : Loo(H) — > Loo(H)- 

From now on we will mainly restrict our discussion to the planar 
case (though we will briefly mention other cases, too). 

As said, our RWwIS is to mimic Lorentz disk process in R 2 . Since 
in the two particle process the energy, i.e. the rate of the particle will 
also change, it is appropriate to include this rate among the internal 
states. Concretely, we will have H — S x X for the set of internal 
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states which now also includes the rate A. Here S = R/Z stands for 
the direction u — of the velocity of a particle and X for its speed 

A = \v\. (Here A G I = [a,b], < a < b < oo; A will be, of course, 
conserved in the absence of interaction). 

Without interaction, the generator for a single random walker can 
be decomposed 

Q = Q s <g> id x 

Also P x {{u,k),.) = Px{u,.) <g> 5\. Thus Q s = E xe z d \{o} p * is indeed 
the stochastic kernel on S. 

2.1.2. Interaction: the collision operator. 

Next we define the collision interaction. Let £\ = (rj\,e\), i = 1,2 
be two RWwIS. 

Whenever n] ^ tj^, the joint generator of the two Markov pro- 
cesses is the product of the two individual generators modeling two 
independent Lorentz processes. On the other hand, when n] = 
nf (= x), then 

P(# + = (x + zW+),8+ =(x + z 2 ,v 2 + y,v\ e A\v\ e A 2 

\Zl = (x,vl),g_ = (x,v 2 _)) 
= C zlrZ2 (v 1 _,v 2 _,A\A 2 ) 

is the collision kernel. We assume that C satisfies conservation of en- 
ergy: (v^_) 2 + (v 2 ^) 2 = {v\) 2 + (v\) 2 (momentum is not conserved 
since the collision kernel contains averaging over normal of impact, 
see below). Thus 

C^CpL, &?.,.,.) = C z i z 2(AL,mL, «?_,.,.,.) 

where Al is the precollisional speed parameter of the first random 
walker (that of the other one is determined by energy conservation). 
For convenience, we will always use the speed of the first walker to 
describe the energy partition between the two particles. 

We can and do assume that (v 1 ) 2 + (v 2 ) 2 = 1. Therefore the state 
space of the two particle process is isomorphic to (Z 2 x S) xl 
where I = [0, 1] (in what follows we always assume 2E = 1). It 
is worth noting that the concrete form of the collision kernel for the 
mesoscopic version of two disk model is calculated in Appendix A 
of IIGG 0811 . 

Warning: the pair £ 2 ) is not a RWwIS on Z d x Z d anymore since 



6 



ZSOLT PAJOR-GYULAI, DOMOKOS SZASZ 



translation invariance is hurt on the diagonal. 



2.1.3. Molecular chaos. 

(1) If we recall that our model is to mimic the two disk process, 
we note that the deterministic law driving the collision does 
not only involve v l _ and but also an angle describing the 
positions of the two disks relative to each other. As frequently 
in the literature, we assume that the distribution of this angle 
is uniform, and averaging over it gives our stochastic collision 
operator defined above. 

(2) As we will see later (cf . Theorem|3]), short inter-collision times 
are extremely rare asymptotically, so the joint law of the di- 
rections of incoming velocities will approach an equilibrium 
distribution. The particular form of this law is not important, 
the point is that by averaging over u x _,u 1 _, we will use the 
mesoscopic collision kernel 

C z i z %{X}_ -> k\,u\,u\) 

This shorthand notation means that we will only use Al from 
the precollisional data to compute the postcollisional veloci- 
ties. 

2.1.4. Summary of the model. Finally, we have as our object of inves- 
tigation 

/t = (fl / g / A ( ) = (^4rf/ef/A t ) 

Introduce the notation A = A + yjl — A 2 , i.e. the sum of the two 
rates and let A, = A if i = 1 and A, = yl — A 2 if i = 2. Then the 
dynamics can be summarized as the following. First, 

&Ut+dt = }t\A t = A) = 1 - Xdt + o(dt) 
For x\ +dt 7^ x\ and x\ ^ xj, 

n?t+dt € {x] +dt } x A,^ dt = f t -\A t+dt =A ( |e| = (x\,u]),A t = A) = 

= AjPj _,(u\,A)dt + o(dt) 

However, when tj] = t]j = x, 

P(g +dt = {x + z 1 } x A',A t+dt G A 3 |d = (x,u\),A t = A) = 

= XC z i z2 (\,A 3 ,A 1 ,A 2 )dt + o(dt) 

Everything else is just o(dt). We launch the process from the initial 
state Jo = (0,u},0,Mo,Ao). 
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Denote the time of the first jump after t with 
t^t)=m£{s>t:J B ^Jt} 
Our main result will concern the four-tuple 

on the set (Z 2 x S) 2 . We will prove weak convergence on the space 
(K 2 x S) 2 endowed with the metric 

2 

d^x 1 , u 1 , x 2 , u 2 ) , (y 1 ,1V 1 ,y 2 ,w 2 )) = £^(|x ! — y l \ + dsiu^w 1 )) 

i=l 

where dsiu^w 1 ) is the length of the shorter arc joining u 1 and w l on 
S. 

For convenience, we define the arithmetic operations on (R 2 x S) 2 
by 

(x 1 , w 1 , x 2 , u 2 ) + (y 1 , w 1 , y 2 , w 2 ) = (x 1 + y 1 , w 1 , x 2 + y 2 , w 2 ) 

c(x 1 ,m 1 ,x 2 / m 2 ) = (cx 1 , w 1 ,cx 2 , u 2 ) ceR 

2.2. Main result. 

2.2.1. Notations and conditions. 

We start with arbitrary dimension d > 1. For a single RWwIS, in- 
troduce the operator valued expected displacement and further the 
analog for the covariance 

(2) Mi = xiP x Zi, m = xix m P x \<l,m<d 

where x\ = {x,e{) and from now on P x = P|. It is easy to see that 
e.g. 

(M / l)( M )=E((^ (t) ) / -(/ /t )zh=«) 

where 1 G Loo ([«,&]) is the constant 1 function on [«,&] and ffj(f) is 
the one walk variant of the above ffj (it will always be clear from 
the context which is to use). Higher conditional moments can be 
defined analogously. Due to the bounded range condition below, all 
these moments are finite. 
Conditions on the RWwIS 
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(1) (Spectral gap) For the operator Q$ : Loo(S) — > Lco(S) 1 is a single 
eigenvalue and the remaining part of its spectrum lies in a circle 
\z\ < 5 < 1. Consequently, for the operator Qg acting over Ai(S), 
the space of finite measures over S, 1 is also a single eigenvalue 
whose eigenfunction is the unique stationary probability measure 
p, i. e. QgjO = p (this can also be written pQ$ = p). 

(2) (No drift) For VI < I < d 

(p,Mil) =0; 

(3) (Bounded range) P x = if\x\ > 1 (for simplicity); 

(4) (Nonsingularity of the asymptotic covariance matrix) Along the 
lines of the LLT of IK Sz 8311 , for the diffusive limit of the displace- 
ments of the RWwIS, the asymptotic covariances of the coordinate 
vectors are equal to 

a lm = (p,Z /w l) - {p,Mi{Qs ~ /) -1 M M l) - 

-(p,M m (Q s - -0 -1 l) 

(1 < l,m < d). It is assumed that the matrix {cim)\<i m <d zs 
positive definite. 

Conditions on the collision kernel 

(1) (Ergodicity) For A C X set g(\-,A) = Ez 1 ^ 2 Qi 2 2 (A, A, S,S). 
We assume that the homogeneous Markov chain defined by this ker- 
nel is an ergodic Harris chain (cf. MR 841 ) with stationary distribu- 
tion distribution p s . 

(2) (Bounded range) Also C z i 2 2 = 0/or \z l \, \ z 2 \ > 1 

2.2.2. Main result. The main result of this paper concerns the limit 
distribution of the two interacting random walkers described in sub- 
section 2.1. 

Theorem 1. For every initial distribution of (£g, £g), the density function 
of the weak limit law of J/ \Tt in ((1R 2 x S) 2 , d) exists and is equal to 

where p is the stationary density of the internal states on S. 

2.2.3. Higher dimensions. In higher dimensions, the treatment is es- 
sentially the same as well as the result with one exception. If F is the 
probability that two particles starting from the same place will meet 
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again, then instead of p s we have to use the mixture 

oo 

p Xo (A) = (l-F)^g n (X 0/ A)F n 

where g n is the n step kernel (Note that F < 1 if d > 3). 

Note that F depends on £q and £g but for the same reasons as 
above, we can assume that it is in principle possible to average over 
them with respect to some certain distribution. This problem is strongly 
related to the question when one asks what the distribution of the in- 
ternal states is at the first return to the origin in the case of a single 
random walker. We do not discuss this in further detail. 

2.3. Methods. 

2.3.1. Local limit theorem for RWwIS's. Our first step will be to gener- 
alize the local limit theorem described in [KSz 83] to continuous time 
and continuous internal states. Since now we only investigate the 
collision-free motion of one particle, the velocity magnitude A will 
be constant, thus we will consider it as a parameter and the internal 
state space will be S. 

Theorem 2 (Local Limit Theorem). With the assumptions in \2.2.1[ namely 

(i) Qs is ergodic and aperiodic with stationary distribution p; 

(ii) (|0,M Z 1) = Ofor every 1 < I < d; 

(iii) The matrix a = {v\m)\<i im <d whose elements are 

trim = (p,tr tm t) - (p,Mi(Q s - l)- l M m t)- 

-{p,M m {Q s -I)- l t) 



is positive definite, 
we have for every ACS 

h tiX (A\£ = (0,u )) 



p(A) ( x 



O ((At) 



-(d+l)/2 



where 

h tiX {A) = P(& e {%} x A) 
and g a denotes the density function of the d dimensional normal distribu- 
tion with mean and covariance matrix o~. The remainder term is uniform 
in uq. 

Remark 1. For the sake of generality we will prove this theorem without 
the bounded range condition only assuming that all moments like © are 
finite and the minimal lattice C determined by the jumps of the RWwIS 
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coincides with Z d (cf. local CLT of |KSz 8 3]). (This condition is called the 
triviality of arithmetics.) 

We will show that - for our case d = 2 - this implies 

Corollary 1. Assume d = 2. If r = inf{f > ffj(0) 

(3) 

F A , M0 (0 = P(r < tfo = (0, mo)) = 1 - + 

ztfere f/ze remainder term is uniform in uq. 

Recall the definition of slowly varying functions 

Definition 3. A positive function L(t) defined on R + is slowly varying at 
infinity if 

^Ut v ceR+ 

Clearly the tail of the distribution function F\ /Uo (t) is a slowly 
varying function, which proves to be a crucial property later. 

By Definition [2j A is just the speed of the random walker and the 
property 

(4) F A;Mo (f) =F 1/Mo (At) 

is evident by rescaling. 

Remark 2. As it will be shown in Section \2.3.2\ these excursions are very 
long, so it is plausible to assume by Theorem\2]that instead of ©, we only 
have to deal with the family of functions 



f a(0 = J FA,u (t)dp(u ) 



The existence of this average is granted by the uniformity of the above 
expansion of F\ /UQ in uq. Rigorously, for every uq, (1 — F\{i))/(1 — 
F\,u Q (t)) — > 1 uniformly as t — > co, and this is sufficient in the sequel 
(cf. the proofs in the appendix). 

Corollary 2 (Central Limit Theorem). With the assumptions of the local 
theorem, 

[%> e tj N d (<d,\o-) x p in (K d x S,do) 

where J\f d (0, Acr) is the d-dimensional normal distribution with mean 
and covariance matrix \o~. The metric do is the sum of the euclidean metric 
on M. d and the previously defined d$. 



ENERGY TRANSFER AND JOINT DIFFUSION 



11 



2.3.2. Scaled Type Renewal Processes. Since the speed of each individ- 
ual particle is conserved between consecutive collisions and the as- 
ymptotic inter-collision times satisfy @, the idea naturally arises 
that one should deal with a renewal process, where the renewal 
times come from a one-parameter family of distributions like (5). 

Consider a family of distribution functions F\ : A € [a, b] C IR+ 
with positive support and assume that VA G [a,b] F\(t) = F(Af) for 
a nondegenerate distribution function F. The corresponding random 
variables are denoted by X A . 

Definition 4. Suppose (Aq, Ai, A2, . . . ) € [fl, fc] N . TTzen the sequence 
S n = E" =0 Xa ; . : n = 0, 1,2, ... is called a scaled-type renewal pro- 
cess (STRP) if X Aq , X Al , X Az , ... is an independent sequence of random 
variables such that V/ £ N the distribution ofX\. is F\.. 

STRP with slow tail return times 

As it was previously mentioned, in the model under investiga- 
tion the return times are very long, more specifically they satisfy the 
slowly varying tail property. 

When Vz A, = 1 (or any constant), there are many well known 
results, among which the one revealing the core of the phenomena is 
the following (cf. IIHM 91IH . 

Theorem 3. Let X, > be random variables with common distribution 
function F. For every k > 0, the following three statements are equivalent. 

• 1 — F(x) = L(x), where L is a slowly varying function 

p 

• X n _ k _ 1/n /X n _ k/n -)• as n 00 

• Li=iXi,n/X n -k,n -A 1 as n -> 00 

where Xj iTl is the ordered statistics from Xj i = 1, .., n. 

In other words, the largest return times dominate the whole process. 
In fact, we do not need this theorem neither its generalization (how- 
ever intuitively it should hold), the important result is related to the 
age and residual age process. 

Now consider the STRP with an arbitrary sequence of parameters 
A = (Aq,Ai, ...) from [a,b]. Let Nt = max{n : S n < t}. The random 
variables 

Yf — t— S;v f Zf = Sw t+ i — t 

are called the age and the residual age respectively. For these quantities, 
we have 



12 



ZSOLT PAJOR-GYULAI, DOMOKOS SZASZ 



Theorem 4. IfF has a slowly varying upper tail, then as t — >■ oo, 



This theorem was proved in a more general setting in BP-GySz 2010 [ 



but in our case a much simpler proof is available which we present 
in the appendix. 

Markovian Renewal Processes 

Definition 5. A scaled-type renewal process S\ , n = YJj=o X\. : n = 
0,1,2, .. . is called a scaled type Markovian renewal process (STMRP) 
if\o, Ai, A2, ••• is a homogeneous Markov chain with values in [a, b] and 
for every realization X\, A2, • • • of this Markov chain S n = YJj=o : n = 
0, 1, 2, . . . is a scaled-type renewal process. 

The essential properties of such processes are well described in the 
literature when the waiting times have finite means (cf . the introduc- 



tion of |P-GySz 20101 for further reference). However, it is clear from 
[I] that we are now facing the infinite mean case which seemed un- 



touched before the authors established results in [P-GySz 2010 [ with 
further restrictions on F. This was done in a more general setting, 
here we only present what is necessary for our current purposes. In 
our case when F is slowly varying, the more complicated machinery 
of the cited paper is not necessary and just as in the case of Theorem 
IH we give a much simpler proof in the appendix. 

Let g(A_, A+) be the transition kernel of the Markov-chain. Sup- 
pose that this is a recurrent Harris chain with stationary measure p s . 

The expectation of X\ is denoted by U\ = u/ A. whenever u = 
f °° xdF is finite. We repeat that the parameter interval is chosen so 
that < a < b < 00. 

Let N t) \ Q denote the number of the renewals occurred before time 
t (including the one at t = 0) with initial parameter value Aq, i.e. 

(5) N tM = inf{n : S Aq , m > t} 

and let U\Jt) be its expectation. 

A classical question is: what is the "type" of the current renewal 
at time t, i.e. what is the distribution of the parameter A. Denote the 
corresponding measure conditioned on the initial parameter value 
A , by <E> Uo , i.e. 

<S> tM (A) = F(A Nuo _ 1 eAC[a,b]) 
By investigating the asymptotics we get 
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Theorem 5. Ifl — Fis slowly varying, then 

lrmO ao (A) =ps(A) 

f— S-oo 

Intuitively, this means that the waiting times are so similar in a 
probabilistic sense, that the rescaling does not matter asymptotically. 
This means that the process behaves analogously as if the waiting 
times were iid. 

2.3.3. Sketch of proof of the main result. For the convenience of the 
reader, here we present the main ideas used in the proof of Theo- 
rem m 

(1) First note, that the spatial difference of the two random walk- 
ers $ = (n} M - n] M , e), e], A t ) is again a RWwIS (modulo 
the origin, where the collision kernel also comes into play) on 
the state space 

Z 2 x S 2 x I 

where I is as in Section 2.1.2. The collision of the two particles 
corresponds to the return of this walk to the origin and we 
can use our results developed for the return times since the 
behavior of the first return is not effected by the dynamics at 
the origin. 

(2) At time t, the relevant information for our goal is the state of 
the process at the last collision before t and its history since 
then. 

(3) It can be shown that, by dividing by \fi, the location of the 
last collision before time t goes to zero in probability. Due to 
Theorem HI the amount of time elapsed since the last collision 
dominates the whole process. Consequently, we only have to 
treat two RWwIS which evolve conditioned on not meeting. 

(4) Now it is clear that the limit distribution will be a mixture 
according to the value of the outgoing A at the last collision 
before t. Due to recurrence, the particles will meet infinitely 
many times and therefore the asymptotic distribution will be 
p s by Theorem |5j 

(5) Finally, we have to derive the limit distribution of two inde- 
pendent RWwIS conditioned on not meeting. The uncondi- 
tional limit is the product of the independent limits deter- 
mined by Corollary |2] and we will show that the conditioned 
one is the same since the condition becomes irrelevant asymp- 
totically. (This is the point where our method brakes down in 
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one dimension). This can be shown by defining an appropri- 
ate random time which behaves like stopping time (although 
it is not), so the strong Markov-property can be used (Lemma 
[5}. Since this quasi-stopping time is very small in a certain 
sense, we can obtain the desired result. However, we will 
choose a different approach and refer to the correspondig re- 
sult for ordinary random walks. 

3. Proof of Local Limit Theorem and of related results 

We will follow the main ideas outlined in MKSz 8311 but for the sake 
of self containedness, we conduct the whole proof except for some 
tedious calculation. 

The transition operator for the discrete time RWwIS is 

T : Loo(Z d xS)4 Loo(Z d x S) 

with 

(Tf)(x)= £ Pyf(x-y) 
yez d \{0} 

for / € Loo(Z d x S). Clearly if T ; is the time when the z'th transition 
occurs i.e. 

T = T t = inf{f > 7f_ a : f]t + 

then 

h T I+1 A-)= E P ?T uX -y(.) 
yeZ d \{0} 

Note that the dual is M. (Z d x S) i.e. the signed measures of bounded 
total variation on Z d x S. 

As usually in case of limit theorems, we use spatial Fourier trans- 
forms 

/: [-71,71)*^ MS) f(s)= £ e^f(x) 

xez d 

For the operator T, we have 

(Tf)(s)= £ &*\Tf){x)= E E e^P y f(x-y) 
xez d xez d yez d \{0} 

which by the change of variables y' = y and x' = x — y further equals 

E E e x ' (s ' z ' +3/) V(*0 = E elW) fy E e i[s ' xl) f{x') 
y'ez d \{0} x'ez d y'ez d \{0} x'ez d 
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that is just oc(s)f (s) where we introduced the operator valued Fourier- 
transform 

cc(s) = £ *W)p v 

y'6Z rf \{0} 

By induction, it follows that 

(6) (IV)(s) = «»(s)/(s) 

Returning to the continuous time, note that by Definition |2j 

{Ti-T^}^ 

is an i.i.d. sequence with distribution EXP(A), thus the number of 
transitions occurred up until time t has distribution POI(At). Using 
this, we have for ACS 



CO 



x,A 



1 ft r n , s 00 (At\ n - — ~-— 



(7) h, x {A\h = (0,«o)) = <T A * £ ^((T*)"(^ Mo )) 

n=0 n - 

Using the Fourier-inversion formula, this equals 

1 / / „-/l5,v) £^ ---^ ' ' 

n=0 

which after using ((6]) becomes 

1 / / .,-/{. -vV) £^ 

n=0 

\nYJ-n"J-n U ° 



(2n) d 



J-n J-tl „_n m 

1 



since ^ ^ = <W 

Proof of Theorem\2\for d — 1. In addition to © we introduce the third 
moment 

s = £ x 3 p x 

xez d \{0} 

following UN 091 . Using these moments 

2 ■ 3 

(8) a(s) = Q S + isM — ^-L — -^-S + o(s 3 ) 

Using a straightforward generalization of Theorem 2.9 in Chapter 
VIII in UK 6611 we have that its largest eigenvalue x( s ) has a similar 
expansion 

X(s) = l + r 1 s + jS 2 + r -^s 3 + o(s 3 ) 
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The coefficients are 

n = r 2 = -(p,H)+2(p,M(Qs - I)- 1 Ml) = -a 1 
as calculated in [K Sz 8311 and 

r 3 = i (3(p # Z(Q s - J^Ml) + 3(p,M(Q s - /^Ll) - (p,El)) 
was computed in IIN 091 . The largest eigenvalue of k(s) = e a *( s ) _1 is 

e^)" 1 = 1 + ns + (r{ + r j)s 2 +(rl+ r -^ + |) s 3 + o(s 3 ) 

by elementary calculation which after plugging the above expres- 
sions becomes 

(9) l-^ + | S 3 + o(s 3 ) 

We mention that if there would be drift, one should use the above 
formula with r x = i{p,Ml) (cf. IKSz 831 ). 

Let <p(s) be the eigenvector corresponding to the eigenvalue x( s ) 
with (p,cp(s)) = 1. Introduce the operator 

P v f = (pj)cp(s) /GLoo(S) 

It is easy to see that for g G .M (S) 

Since every moment is finite, the perturbation is analytic and e a ( s )~ 
is continuous in s. Consider the operator valued function 

R(s)=e^- 1 -e^- 1 P (p{s) 

which is again continuous in s. Since by assumption there is a spec- 
tral gap for Q s , 

\\R(0)\ \ = ||e Qs_1 -Pi|| < 1 

Continuity implies | \R(s) 1 1 < 1 for sufficiently small s. 

Since the eigenspaces depend continuously on the perturbation at 

s = 0, 

e^ 1 =(l-^f + T js 3 + o(s 3 )) (P t + o(l)) + R(s) 
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and thus e 



At 



(7 Z S Z f3 
1 - -^r- + — ~ — ,„ + 



can be decomposed 

„3 N N Af 



2Ai 6 (Af) 3 / 2 V(Ai) 3/2 



(P! + o(l)) + 



(10) 



Af 



where the last term is exponentially converging to zero. Now we 
will show that 2n^/Xt times 



2tt 



71 



(11) 



g(A) 

\/2nMc 



e 2\ta- 2 ( 1 



z>3 x(3cr 2 Ai — x z 
T cr6(Af) 2 



goes to zero. By x/(s) = —isf'{s), after some elementary calcula- 
tions, one can get 



6 u 6 (A() 2 

r oo 



V2nJ-co V 6 / 

plugging this back and making the change of variables s' = a/aTs we 
have 



nsJJd 

e 

-nVXi 



p(A)j_ 



— 1Y— —as' 



2J2 



r 3 s 



13 



6 y/Xt 



ds' 



By the triangle equality this integral is bounded from above by the 
sum of the following four terms (we drop the prime) 



h = 



s|<(Af) £ 



Wis)- 1 )*, 

h = p{S)C f 



p{A)e 



1 + 



T3 S° 

6 x/AT 



ds 



'ds 



h 
h 



s\>(MY 



L 



(At) e <|s|<7VAt 
7V / Ai<|s|<7rV / Ai 



ds 
ds 
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for < e < 1/6, C = J (l + and Af > 1. Again by elementary 
calculations, 



a 2 s 2 



+ 



''3 



2At ' 6 (Af) 3/2 
= e 2 



+ o 



(AO 



3/2 



Af 



1 , r 3 S a 
1 + — -= +0 

6 VAi 

From (jlO|) now it is clear that Ji = o(l/vAi). Using the usual 
upper bound for the tail of the Gaussian function, it is easy to see 
that I2 = o(l / vAf) as well while the trivial arithmetic condition - 
similarly as in [KSz 83J - implies the exponential decay of I4 with 
growing At. If 7 is chosen small enough, then from ([9]), it follows 
that 



< e 



With this 

h =VTt 

< 



(Af) e ~ 1/2 <|s|<7 
(Af) e - 1/2 



»o 



\s\ < 7 
ds < 



(7 2 S 2 A< 



4 ds 



< s <7 



4 ds = 



'(Af) e <|s|<7VAt 

as before, so Ij + I2 + J3 + 4 = J an d the proof is ready □ 

Proof for d > 2. The multidimensional case is a straightforward gen- 
eralization. For expansion of the largest eigenvalue 



M-i 



[s,crs 



d d 



r 3,i,j,k 



2 +EEE^f Zs ' s J s fc + °(i s i 



i=ij=ik=i 

while one has to prove the convergence of 

(Af) d/2 



{2n) d 
P( A ) 



r ... r e -M e ^(s)-vs uo ds 

J — 71 J — 71 



00 At(s,o-s) 

e 2 



-z(x,s) 



1 + Xtf(s))ds 



(Af) d/2 

to zero where /(s) is the above term containing the r^j^-s. It turns 
out that the term containing f(s) is just O (l/ y/Xtj and Theorem|2] 



follows. 



□ 
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Proof of Corollary^ We will obtain the formula © through the dis- 
crete time process. The discrete process was defined in Definition [1] 
As a straightforward generalization of Theorem 2 of IN 091 using the 
elements of the previous proof, we have 

KM) = F (£« e {x} x a|£ = (o,m )) = 



For d = 2, the remainder term is 0{n 3 ^ 2 ) which means it is sum- 
mable. From the proof of Theorem 6 of IN 0911 , we know that 



/„ = P(t > n\rj Q = 0) 



| I log log n 

log n 



log n 



with the remainder term being uniform in the initial state. From now 
C := 2n^/]a\. By the same argument that led to ©, 



First, introduce a cutoff 



-At f (A0 M f 



n=0 



L|Afj 



-Af 



(AO" 



< P( |N t -Ai| > ^ ] < 4 



using Chebyshev's inequality. For the remaining n's it can be easily 
obtained that 

MA»)_ 1 + C / i 



logn \.log(Af) 

Using this and the Chebyshev inequality again after some calcula- 
tion we obtain 

y 1 f e -*(*£ = _^_ + o ( 1 "l + 

_Vl n »' log(At) +U lAtlog(At)J + 

+ ( — ^ ) + O f \ I + RT 

The remainder term can be estimated similarly. With some ele- 
mentary calculation again, 



loe loe n 
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and by short computation we obtain 

RT = /loglogW 
I log 2 (At) 

Finally one can observe that this last one is the slowest error term 

so 

i-F(t)=fffi+of''*y t > 

' log(Af) ^ l g 2 (Af) 
as desired. □ 

Proof of Corollary^ The result could be derived directly from the lo- 
cal theorem but in one dimension, it is simpler to let t — > oo in ((T0|) to 
obtain 



AfU*( -4= P _<d£ 



where by definition P|£ M = p. 

The multidimensional theorem can be obtained similarly □ 



4. Proof of Theorem 1 

First we want to transfer our results for the returns of one random 
walker to the origin to the collisions of the two particle case. Note 
that the spatial difference of the two random walkers form a one 
particle RWwIS £j = {rj} XQ — e\, ef f A;) on the state space 

(Z x S) 2 x I 

except for the origin where the collision kernel spoils the translation 
invariance. The component in I is the speed-parameter Af of the first 
particle, while the other one's is determined by the fixed total energy 
(2£ = 1). However the speed-parameter of this difference-process is 
the sum of these two individual parameters (cf. superposition of 
Poisson processes) 

A t+v /l-A 2 

This is clearly bounded away from zero except for the trivial zero 
energy case, so we can use the results in Section 12.3.21 

Although the behavior is different in the origin, the tail of the times 
between consecutive visits to the origin of £| (collisions) will still be 
slowly varying. 

Thus if r(t) is now the time of the last collision before t, i.e. 

T(f) = sup{s < f |3e > : r\\ M £ r]\ M , r\\ M = y]\ m s-e<u<s} 



ENERGY TRANSFER AND JOINT DIFFUSION 21 



then r(t)/t 4- Oby Theorem H 



Now we have the following decomposition of the process 

= h-r(t) , Jr(t) 

Jo V * V^W) 



(i 2 ) Jl 



+ 

r(f) 



Jo 



where we indicated the starting states in subscripts. The last term 
can be dealt with using Theorem 4.1 in [B 68J and the following 

Lemma 3. For the common place of the last collision rf T ^ = Vr(ty 

Proof. Pick e,S > and let f be large enough such that P(r(i) > 
et) < <y. Then 



P 



Vt 



> e) < J F(\n c ut \ > eVt\r(t) = ut)u T ^(tdu) +S 



where }i r (t) ^ s * ne measure generated by r(f). 
Since 

{|< t | > eVt} n {T(f) = ut} C {|^| > ev 7 *} n {T(f) = ut} 
the integrand can be bounded from above by 
(13) F(\nl t \ > eVt\r(t) = ut) < 

< supPd^fl > eV~t\nl t = t]l t ,A ut = \,e l ut = u l ) 

A,m' 

The inequality is due to the fact that 

{r(t) =s} = {rjl = n 2 s } n {irf + n\ ut (s,t]} 

where the second event only affects the distribution of n\ through 
the internal states. Obviously the condition n\ t = n\ t is restrictive 
spatially so we can further estimate (fi~3)) by 

supP(|^f| > eVt\A ut = Ke l ut = u 1 ) < 

\,u' 

1 

< ^-supD 2 (|^||A Mt = X,e l ut = u 1 ) 

€ 1 A,u' 

using Chebysev's inequality. Note that by the diffusive nature of the 
process, the second moment is monotonously increasing with time, 
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and it is the largest if r\\ posesses all the energy throughout the whole 
process up to time ut. Thus, the further bound can be obtained: 

isupD 2 (|^||A s = 2E s e [0,6t],4t = u ') = 

u' 

= 4tSupD 2 (|^||A s = 2E sg [0, St], e} t = u 1 ) 
e t «i 

since the second particle is then standing still at the origin. This vari- 
ance is nothing else but the variance of a single RWwIS rj$ t with rate 
2E and its internal state conditioned to be u 1 at St. 

By the bounded range condition, the above variance can be esti- 
mated from above by the variance of an ordinary continuous time 
random walk, in which the one step variance is 1. It is well known 
that the variance of such a process at St is 2EStcr for some constant &. 

Using this and that u T ^([0,tS]) < 1 we have that 

P {^j^- >^j< S(e~ 2 2Ea + 1) 

for some constant K independent of S. Since S is arbitrary, the proof 
is finished. □ 



As the next step, note that y/ (t — r(t))/t A 1 by Theorem H] and 
by Theorem 4.4 in [B 68] we only have to prove the weak limit of the 



remaining term in (|T2|) . 



This limit is nothing else but the joint limit of two continuous time 
RWwIS starting from the origin and not meeting once they depart. 
The distribution of the energy between these two are according to 
the distribution of A T m and by Theorem|5l we have 

P(A T(f) e AC I)^p s (A) 

Thus we are ready if we can show that 

Lemma 4. In d = 2, the joint law of two independent RWwIS (starting 
from the origin) conditioned on not meeting once they depart is the product 
of the independent one particle limit-laws, i.e. for continuity sets (of the 
appropriate measures) A\ and Ai, 

nl + nl se[t ( (o),i] 

(is) ^^ F (UA €Al 
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Before we proceed with the proof, we establish an identity first. 
Introduce 

T(t) = inf{s > : n] = n 2 s , V l / rj 2 u t^s) < u < s + t} 

Note that if there is a constant period of length t at the same site be- 
fore a "depart and not return", this definition of T(t) gives its start- 
ing time. Note that P(T(f) < oo) = 1 and that (T(f) > a} C {t(a) < 
i} and thus 

,(M > .), P (^ < I).o 

In other words, T(f)/f 4 0. 

Let Jf d denote the non-interacting two-particle system which makes 
it's first jump according to the collision kernel and note that then 

[Jtlvl^Vs se[t n {0) f t]]^\ff\r,l^r,l s G [£fj(0), f]] 
We have 

Lemma 5. For any A £ £>((Z 2 x S) 2 ) 
P(i? e AI^V /y 2 se[f fj (0),f]) = 

= (i - e - Af )F(/? (t)+t - Jrw g A l4 (f) = 4) + 

A = A + Vl - A 2 . 
Proof. Assume first that ffj (0) < t and set 

L s = Vue[o*-t]{vl = nl>nl £ y 2 >tf)W <r<u + t} 

and let \it denote the distribution of T(t). Then 

F (/r(0+t ~~ Jr(f) G A l e T(t) = e o) = 

= / nilit - ff e A\T(t) = s,ei = 4)df*t(s) 

which equals 

/ P(/i d +f- Ts e A|Lg, //s = j^i ^ J/, 2 rG [f^sLs + f],^ = 4)d^ t (s) 

Since our process is Markov, L c s is superfluous while the remaining 
integrand is time-translational invariant, so by P(T(f) < oo) = 1, 
rjl = rjQ = and the way we defined the substraction, we obtain 

nft d e A\rj l r £ rel^f]) 
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On the other hand if ffj(O) > f, then T(t) = thus we get 



P(/l d G A|ffj(0) >t) = S 



Jo 



□ 



Proof of Lemma^ By assumption for every P\ -continuity set A C 

R 2 xS 

id 




e A 



1,2 



(Recall that A, = A for i = 1 and A, = vl — A 2 for z' = 2) where the 
P\.-s are determined by Corollary [2} thus all we have to show is 



(16) 



P 



' fid 
J T(t)+t 



fid 

h(t) 



-'T(t) 



Although this could be veryfied directly, we choose a different ap- 
proach. For ordinary random walks, the authors - generalizing a re- 
sult of Bolthausen - established the desired result in a functional con- 
text (cf. Corollary 1 in and remark (4) in |P-GySz 2010b [). For eco- 
nomicity, we omit the proof of an invariance theorem of the RWwIS 
(although it is not by any means harder than the invariance princi- 
ple for the ordinary RW) and the obvious generalization of the cited 
result. □ 



5. Remarks 



(1) The aformentioned direct proof of Lemma|4]is based on T(t) / 1 
0. It suffices to show that 



P 



TJ 



T(t)+t *7r(t) 'it 



Vt 



> e 



-T(t) 



-^0 1 = 1,2 



Similarly as in Lemma|3l 7^L^ / \/t — >, so we can drop it from 

the above formula using the triangle inequality. Then, we 
can argue that T(t) being small implies that the difference of 
\f]r(t)+t ~ Vt \ is small enough such that it converges to zero in 
probability in the scaling limit. 
(2) For ordinary random walk there is a variant of the local theo- 
rem which is a better spatial estimate (cf. P7.10 in [S 76J). The 
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corresponding theorem for RWwIS is 



xez d 



E I' 



2 ht, x (A\£ = (0,u )) 




This gives the limit of the mean and the variance of the abso- 
lute value. 

(3) For treating the deterministic model, the realistic alternative 
is to rely upon the averaging method of HChD 091 . Indeed, 
between two collisions of the disks there typically occur long 
collision sequences of the particular disks with the periodic 
configuration of fixed scatterers. During these long intervals, 
their orbits become approximately Brownian and their veloc- 
ities and the normal of impact incoming into a particular col- 
lision of the two disks correspond to an equilibrium distri- 
bution and finally their outgoing velocities from the collision 
can be calculated analogously to the collision operator ap- 
pearing in the derivation of Boltzmann's equation for a hard 
disk fluid. We plan to return to the deterministic model in the 
future. 
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Appendix A. Proof of STRP and STMRP results 

A.l. Preparatory facts. To prove our results we need the so called 
Abelian-Tauberian theorems (see ||F 70] XIII.5). Before presenting 
them, introduce the Laplace transform of F: 

rCO POO 

m( z ) = / e~ zx dF(x) = z / e~ zx F(x)dx z > 
Jo Jo 

where the second equality can be obtained by partial integration. 
Clearly by the scaling relation, <p\(z) = <p(z/X). 

Fact 1 (Feller). Let Hbe a measure on R + , k(z) = J e~ zx dH the Laplace 
transform wrt it and H(x) = H( [0, x] ) / Then for p > 0, 

x(t/x) „ 

K{l/X) 
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and 

— ^ -HP X -> 00 

zrapZy eac/? of/?er, moreover in this case 

(17) k(1/x) - H(x)T(p + l) x^oo 
A popular reformulation of this result is 

Fact 2.I/L is slowly varying in infinity and < p < oo, then 
k(1/x) ~ x p L{x) x — > oo 

and 

1 

H(x) ~ =7 ^x^Lfx) x — > oo 

v ; r(p + i) v ; 

implies each other. 

We also need a technical result (cf. |F 70] VIII.9 Theorem 1) 
Fact 3. IfK(x) is a slowly varying function at infinity and /3 > — 1 then 

J*yPmdy = (t^6 X ^ K ( : 

A.2. MRP results. As mentioned in Section 2.3.2, we will investi- 
gate the asymptotic behavior of <3>^ Ao (A) as t — > oo. By conditioning 
on the first return we can write 

<Z> Uo (A) =l{A eA}(l-fA (0) + 

+ f [ g(A 0/ A 1 )<D i _ S/Al (A)dA 1 rfF Ao (s) 

JO J a 

This is, however, not the usual renewal equation, thus we have to 
generalize the standard results of renewal theory. 

The unique solution among the functions which are bounded on 
bounded intervals is 

(18) <^ Ao (A) = / (l-F A (t-s))U Ao (ds,dA) 



X 



where 



U Ao (f,A)=E^N t/Ao l {ANao _ ieA} ^ 



Note U Ao (t) = U\ Q (t, [a,b]) is the usual renewal function. 
Introduce the measure 

(19) cd Ao (z,A) = [ e- zs dU Ao (s,\) 

J [0,oo] x A 
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Also let S(z, A) be the Laplace transform of <3>f /Ao (A) in the variable 
f. 

Then by ((T8l> , Fubini's theorem, and the product rule of the Laplace 
transform, 

(20) S Ao (z,A) = ! 4> A (z)dtd AQ (z,A) = - f (l-<pi(z/\))(D(z,d\) 

J A Z J A 

V v ' 

where the integration is wrt the measure defined in (|T9|) and 



poo ^ — (0{ — 

(21) <p A (z) = / e- zx (l-F A (x))dx~ 1 

Jo 

Proof of Theorem\5\ We need to show that W Ao (z, A) ~ |0 S (A) as z — )• 
because then the proof is ready by Fact |2] since it implies 

(22) f <P S)Ao (A)ds ~ tp s (A) t^oo 

J 

From this, it is easy to see that <E> f/Ao (A) — > p s {A). 
By the monotonicity of cp, 

(l-cp (£)) d) Ao (z, A) < W Ao (z, A) < (l - cp (£)) d) Ao (z, A) 

We prove that the upper bound is ~ p s (A), then the case of the lower 
bound is trivially the same. 
It is not hard to see that 

d> Ao (z,A) = l {Ao6A } +£(A ,A)<pi (j^-j + 

+ <Pi (j^j g(Ao/ A x )(pi f ^ g(A X/ A)<Ui + ... 
Again by the monotonicity of tp, 

(23) d> Ao (z, A) < l {Ao€A} + £ ^ (-) g"(A , A) 



n=i 



We also know that the convergence of the Markov chain is exponen- 
tial, i.e. 

sup |g M (A , A) - p s (A)\ < C x e-^ n C lr Ci > 

AeB{[a,b]) 
Using this, <(23]) equals 

(24) p s (A) * +KT 

1 - 9 UJ 
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where the remainder term can be estimated 

\RT\ < \t {XoeA} -p s (A)\ + (r) \g n (^o, A) -p s (A)\ < 

^ 2 + c ' i-«-^ (j ) =0(1) "° 



By multiplying (|24|) with 1 — cp{z/a) and using that 1 — cp is a slowly 
varying function as z -4- 0, the proof is finished. □ 

A.3. STMRP with slow tail return times. 

Lemma 6. 

(px(z) ~ L A Qj so also 1 - (z) ~ F A z ->• 

where L A (f) = 1 — F\(t) as t — ^ oo. 
Proo/. Let 

H A ([0,x]) = H A (x) = f 1 - F x {s)ds = (1 + o(l))x(l - F A (x)) + o(x) 

Jo 

by Fact 3, so 

,25) |f^ f 
Due to Factd 

^a(z) t 

and by CL7j>, 

1 - <p A (z) ~ = z | q 1/Z (1 - F A (s))ds z ^ oo 

By (|25|) and the definition of H A/ we can write 

H A (*)~fL A (0 f^oo 

where 

L A (0 = ~ /'(I - F A (s))ds ~ 1 - F A (f) t ^ oo 
t Jo 

which has L A (f) = L(Af) = 1 - F(Af). By FactH 

(26) <£ a ( z )~Il a QJ 1 - <p A (z) ~ L a Qj z^O 

□ 
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Introduce now 

00 

u A (t) = E( n *)f=i f A,(0 

n=0 

where n* denotes the convolution product. 
Lemma 7. 

U A (t){l-F l {t)) -vl t^oo 

Proof. Consider first the case when all A, = A and introduce the func- 
tion LZj(t) = E~=o f l M (0 and its Laplace transform co~ x {z)\ Then 

00 1 1 

<^a(z) = } <P?(z) = r^r ~ 1-7 z — ^ 

z/ 1/z (l-F A ( S ))ds 

Using this with (|26l> , we have 

^ (2) ~II(T7i) 

which is a slowly varying function at zero (which means 1 z) is 
slowly varying at infinity By Fact|2l 

1 1 

(0 ~ t — rr ~ ^ ^~rr * 00 

I*(0 !- f a(0 

so we finally have by 1 — F being slowly varying that 
(27) (l-F(t))U°(f)^l f^oo 

Now we will prove that (|27|) also holds in general. Clearly 

U° a (t) < U A (t) < 11° (F) 

so by (|27)) , we have 

LT A (f)(l -F(f)) -V 1 f-^oo 



□ 



Proof of Theorem® Now pick arbitrarily small e and calculate 

P (*<!-« 

which is 

P (u~ =0 U y6[e/1] {S n = fy} n {X An+1 > f(l - y)} 

oo „1 

= E / (i-fA n+1 (i(i-y)))^((n*)f =1 F Ai (fy)) 
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00 f i 

= E / (i-F(W(i-y)))^((n*)f =1 F Ai (^)) 

n=0 Je 

which has an upper bound by Aj > a > 0: 

f 1 ,w nnn J1T / \ f 1 1 -Fi(flt(l - y)) dU A (ty) 

as f — > oo. Here the first term goes to zero everywhere except y — 1, 
since 1 — F(f) is slowly varying, while we will show that the mea- 
sure wrt we are integrating, converges uniformly to the point mass 
concentrated on zero. To see this, Laplace transform the measure 
dU A (ty)/U A (t) to get 

1 fV z W A (iy) - " A(ZA) - 



U A (t)Jo U A (t) 
due to Theorem [1] which is the Laplace transform of the point mass 
on zero. Therefore, on the domain of interest this measure uniformly 
goes to zero. 

Since e is arbitrary, the desired result follows. The result for the 
residual lifetime can be obtained similarly. □ 
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